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Abstract 

Using a purely hadronic model, we study the charge-symmetry-breaking p-u, 
(T-(5, (T-p and 5-u: mixing amplitudes in isospin asymmetric nuclear matter. 
The basic assumption of the model is that the mixing amplitude is generated 
by nucleon and anti-nucleon loops and hence driven entirely by the difference 
between proton and neutron Fermi momenta and the proton-neutron mass 
difference. We find that the behavior of the mixing amplitude is very com- 
plicated in the spacelike region and quite sensitive to the proton fraction of 
nuclear medium and nuclear density. In particular, in neutron rich nuclei 
(like Pb) and/or neutron stars the mixing amplitudes become about 10 ~ 100 
times as large as those in typical mirror nuclei. 
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Charge independence and charge symmetry (CS) provide powerful tools in organizing and 
describing the multiplet structure of hadrons and nuclei. These symmetries are imperfect 
and, in particular, charge symmetry breaking (CSB) has been studied for a long time [|1|. 
It is trivially broken by the electromagnetic interaction. Other well known sources are the 
proton(p)-neutron(r;,) mass difference and the nucleon-nucleon (NN) potential caused by 
meson mixing, and these sources can be related to the u-d quark mass difference in QCD. 
Thus, the small imperfection of CS provide a unique opportunity to investigate the relation 
between QCD and hadronic and nuclear observables. 

It has long been recognized that the p-uj mixing amplitude would give rise to the signif- 
icant part of CSB in the NN force. Indeed, substantial effects have been observed in the 
e^e' TT+TT^ reaction at ~ (m^ = 782.57 MeV, the uj meson mass) ||l]J^. These 
results allow an extraction of the p-u mixing amplitude, {p\Hmix\^) — —4500 MeV^ [0, 
which seems to be consistent with CSB in the NN scattering lengths the Okamoto- 

Nolen-Schiffer (ONS) anomaly |T],^ and the analyzing powers measured in elastic n-p scat- 
tering 1,1. 

However, about a decade ago, the importance of the momentum dependence of the p-uj 
mixing amplitude was pointed out using several different models It suggests that the 

p-uj mixing amplitude at spacelike momenta is different from its value at the uj pole. A 
QCD sum-rule calculation also gives a strong momentum dependence of the amplitude 0. 
It implies that the NN potential given by the off-shell amplitude of the p-uj mixing is quite 
different from the ones used in the successful phenomenology. 

On the other hand, another analysis argues that the strong off-shell effect is inconsis- 
tent with the observed g^-dependence of p-7* coupling (see also Ref. Moreover, Cohen 



and Miller |]T0[ have emphasized that simple knowledge of the off-shell meson propagator 
is not sufficient to determine the CSB NN potential and that it is important to know the 
vertex function computed from the same theory that supplied the propagator at the same 
time. They have showed that the momentum dependence of the p-u mixing amplitude may 
be absorved into a re-definition of the p-N vertex. However, this problem is still controver- 
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In this Letter, our aim is to study the properties of CSB meson mixing in isospin asym- 
metric nuclear matter. To understand CSB in nuclei (for example, the ONS anomaly, nuclear 
structure of heavy nuclei like Pb, neutron stars etc.), it is vital to know how CSB effects 
depend on nuclear density {ps) and numbers of proton (Z) and neutron {N). Such a study 
has not yet been done in the spacelike region |]T2[. In the present calculation, we use a purely 



hadronic model, i.e.. Quantum Hadrodynamics (QHD), to evaluate various meson mixing 
amplitudes in asymmetric nuclear matter. 



According to Ref. [|I0|, let us introduce a typical scale A^ which separates the short ranged 
effects from the long ranged ones in meson exchange models. Then, the short ranged parts 
are handled by vertex functions, while the long ranged ones are explicitly treated by meson 
propagators. Here A^ is taken to be of order of 1 GeV so that a, p, u and 5 (or oq) mesons are 
explicitly included but heavier mesons are not. (In the present calculation, the pseudoscalar 
mesons are not included because we suppose the relativistic Hartree approximation (RHA).) 

The Lagrangian density is thus written as 

^ = ^QHD-I + i^NNS + i^NNp + SCqT, (l) 
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where Cqhd-i is the usual QHD-I Lagrangian [|T3| which describes the a-uj-^ system. The 
last term in the right hand side is a counterterm for renormalizations, which will be discussed 
below. 

The 5 meson is expressed by the Lagrangian density C]^m5 and it interacts with the 
nucleon through an isovector, scalar coupling 



£int. 
NN5 



(2) 



where 5 is the neutral 5-meson field and gs is the 6-N coupling constant. The nucleon mass 



in matter M* 



p OT n), which is a function of the a and S fields, is now given by 



M(p\ - g^a =F gsS, 



(3) 



with a the a field, g^- its coupling constant and Mj the free proton (Mp = 938.27MeV) or 
neutron (M„ = 939.57MeV) mass. 

The isovector p meson has a vector as well as a tensor coupling to the nucleon. The p-N 
interaction, which is involved in the Lagrangian density Cnnp, is given by 



nint. 
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f _ 
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where and p^ are, respectively, the fields of the nucleon and the neutral member of p 
meson. Here gp and fp are the vector and tensor coupling constants, respectively. 

To study the meson mixing amplitude, we first have to solve the nuclear ground state 
within RHA. In case of the p and uj mesons, since the vector self-energy due to the nucleon 
loop vanishes in RHA ITB,!^], the vacuum polarization does not change the energy density 
of the total system. 

For the scalar mesons, it is necessary to treat the a-5 mixing and the polarization correc- 
tions to the a- and 5-meson propagators simultaneously. In Eq.(|I]), 5Cct involves countert- 
erms such as cr"(5™ (n, m = 1 ~ 4), which remove the divergent pieces of the nucleon-loop 
contributions to the a-d mixing and the scalar- meson propagators. As in Ref. []T^, the 
counterterms are determined so as to renormalize the polarization functions at = and 
M* = Mj. This prescription provides a finite vacuum fiuctuation correction to the energy 
density of the total system. We finally find that in RHA the correction is given by a sum of 
A£p and A^^„, where 



87r2 



M*Mn 



M* ~ 

-^]+Mf{M,-Mr. 



-M^{M,-M:f 



1 'i 2^ 
+ -M,{M,-M:f--{M, 



m: 
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The total energy density in RHA is thus given by 



£tot — S-MFA + A£^p + A£^„ 



(6) 



where Smfa is the energy density in the mean- field approximation (MFA) []T3[. In the present 
calculation, the a, 6 and p masses are, respectively, taken to be m„ = 550.0 MeV, ms = 983.0 
MeV and nip = 769.3 MeV. We determine g^j and g^ (the tu-N coupling constant) so as to 
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fit the saturation condition for symmetric {Z = N) nuclear matter: Stot/pB — M = —15.75 
MeV (M, tlie average of the free proton and neutron masses) at the saturation density 
Pq = 0.17 fm~^. Fixing the ratio of the tensor to vector couphng constants, Cp = fp/gp, 
to be 5.0 [|l^], the vector couphng constant gp is chosen so as to reproduce the measured 
p-u mixing amplitude, {p\Hem\i^) = —4520 MeV^, at the u pole 0. Finally, we determine 
gs so as to produce the empirical symmetry energy a4(= 32.5 MeV) ||16|. We then find 
that g„ = 8.360, gs = 3.257, g^^ = 9.493 and gp = 3.455. This yields the effective proton 
and neutron masses, M* = 679.36 MeV and M* = 680.65 MeV, at po and the nuclear 
incompressibility K = 460 MeV. 

Now we are in a position to calculate meson mixing amplitudes in asymmetric {Z ^ N) 
nuclear matter. There are four types of the CSB meson mixing: p-u, a-6, a-p and 6-uj mixing. 
In particular, the mixing of scalar and vector mesons (the latter two cases) is similar to the 
familiar a-u mixing in matter |ll^ , p!7| . It is a purely density effect and forbidden in vacuum. 
Note that the p-uj mixing yields class III and IV forces while the others provide class III 
forces ||l]]. 

We first calculate the p-u mixing. The polarization function due to the nucleon loop 
has a tensor character and it can be separated into the longitudinal (L) and transverse (T) 
components. They are identical each other in vacumm. The mixing amplitudes for the L 
and T modes are, respectively, given by the polarization functions as 

1 

{p\H^i^\u)L = -|n°° (go, Qs) and {p\Hmix\^^)T = ^[nJi(go, Qs) + ^%{qo, q,)], (7) 
where we choose the direction of the meson momentum as = (go, 0, 0, g^) and 

n;i:;(go,g.) = -^g.gpJ ^Tr[7^G(fc + g)f;r,G(fc)]. (8) 



Here G is the nucleon propagator |[T3| and is a diagonal matrix in isospin space 




As usual, the nucleon propagator can be devided into the Feynman (F) and density- 
dependent (D) parts |T^. 

In the p-N tensor coupling in matter, it is natural to replace the free nucleon mass in 
Eq.(^ with the effective mass M* because it is originally derived by the Gordon decompo- 
sition of the nucleon field in Hence, in Eq.(p|) we use the p-N coupling as 

= 7^ - i4l-t^''^9A and f^t = -^^ + i^a^^qx. (10) 



This replacement of the nucleon mass is equivalent to adding terms such as a^'i/ja^'^di^pp^Tzi/j and 
S^ipa^^ duPp_Tzip (n = 1, 2, • • • , oo) to the original Lagrangian density. However, those terms do not 
modify the a and 5 fields in RHA. 



4 



Because of the isospin matrix t^, the mixing amphtude is given by the difference between 
contributions from the proton and neutron loops. The polarization function which involves 
only the F part of G becomes divergent. However, as in Ref. [0], by virtue of the replacement 
in Eq.(|lOD the difference between proton and neutron contributions becomes finite even in 
matter. The F part of the polarization function is then given by a sum of and nf 
where the former is the polarization function with the vector coupling at the p-N vertex 
while the latter is that with the tensor coupling: 



pu 2 



27r 



dxx{l — x) In 



M*2 
p 



x{l — x)ql 



dxln 



M*2 

n 

Mp*2 - a;(l - x)q^ 
M*2 — a;(l — x)q'f^ 



x{l~x)ql 



pj 



(11) 
(12) 



with C,'^'^ = —g^" + q^if I'^ji- The polarization function involving at least one power of the D 
part of G can be calculated analytically 

In the present calculation, we introduce a monopole-type form factor at each vertex in 
the spacelike region. Thus, we replace the coupling constant as 



9 



9 /{I - ql/A' 



(13) 



1.5 GeV. Here we suppose that A is common to all the vertices. (As 
p!0[| , properly speaking, it is necessary to calculate the CSB vertex from 



with the cutoff A = 
pointed out in Ref. 
the same QHD theory that provided the propagator.) 

In Fig. |T], we present the p-u mixing amplitudes for the L and T modes in the spacelike 
region. We fix go = and vary |g|(= g^). The nuclear density is chosen to be po- Here we 
introduce the proton fraction fp that is defined by Z/ [Z + N). In the figure, in addition to 
the case of symmetric nuclear matter {fp = 0.5), the results for fp = 0.49 and 0.51, which 
respectively correspond to the mirror nuclei of ^^Ca and ^^Sc, are illustrated. The mixing 
amplitude in symmetric matter is very similar to that in vacuum.l However, in asymmetric 
matter the amplitude is much larger than that in symmetric case. The mixing amplitude 
for fp = 0.51 is opposite in sign relative to that for fp = 0.49 because the amplitude is given 
by the difference between proton and neutron contributions. It should be noticed that there 
is a sharp peak around = 2kp {kp, the Fermi momentum of proton or neutron) in the 
mixing amplitude and that the peak structure stems from the D part of the polarization 
function. As an example, in Fig. ^ we illustrate how such a peak is generated in the T mode 
with fp = 0.51. From the figure, we can see that the peak is produced by the cancellation 
of large amplitudes from the proton and neutron loops. In Fig. |], as a function of Pb and 
fp, we show the ratio of the p-u mixing amplitude at the u pole (go = and g^ = 0) in 
asymmetric matter to that in vacuum. Note that in the limit g^ ^ (go 7^ 0) the L mode 
is identical to the T mode. The mixing amplitude at po decreases gradually as fp increases. 
However, at high ps the amplitude changes drastically. In particular, at small fp and high 
Pb the ratio becomes 30 ~ 40. Such an extreme condition may correspond to the case of 
neutron stars 



16,18 



^The present result in vacuum is consistent with the result calculated in Ref. Q. 
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Next we consider the a-S mixing. The mixing amphtude (or the polarization function) 
is given by 



{5\H„nxW) = Il„siqo,qs) = -Wa95 J —--^TT[G{k + q)T^G{k)] 



(2vr) 



(14) 



The amphtude is again given by the difference between contributions from the proton and 
neutron loops. As mentioned above Eq.(^), the a-6 mixing does not vanish in vacuum 



if Mp 7^ Mfj. The normalized F part of the polarization function is given by 11^^ (g^^ 



- / dx [M*' - x{l - x)qt\ In 



*2 



x{l — x)q 



(15) 



The D part of the polarization function can be calculated explicitly. It is similar to the 
D part of the nucleon-loop contribution to the cr-meson propagator, except for the isospin 
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For the scalar-vector (cr-p or 6-uo) meson mixing, the mixing amplitude (or the polariza- 
tion function) is given by 



{v\Hmix\s) = Ilsv{qo,qs) = WsQv^i 



d^k 



Tr[G'(A; + g)f^,G'(fc)], 



(16) 



where s and v stand for the scalar and vector mesons, respectively (if v = u, = 7^). 
As in the a-u mixing [^,0, this amplitude vanishes identically in vacuum and hence is a 
purely density effect. The polarization function can be calculated analytically |1^ . 

In Fig. ^, we present the results of the a-6, a-p and S-u mixing in asymmetric nuclear 
matter. The behavior of the a-6 mixing is very similar to that of the 6-uj mixing. However, 
the a-6 mixing amplitude is opposite in sign compared with the S-uo mixing. By contrast, 
in the a-p mixing the sign of the amplitude changes at ~ 0.63M. This is due to the p-N 
tensor coupling that dominates for large g^- In the top panel of the figure, the a-6 mixing 
amplitude in vacuum is also shown by the dotted curve. It may also contribute to the NN 
scattering lengths in free space ||l|]. 

In Fig. 1^, we illustrate the meson mixing amplitues in Pb [ps = 0.15 fm~^ and fp = 0.4) 
and neutron stars. We suppose that in neutron stars the nuclear density is about 3po and 
0.15 [ll6| , p!8| . From the figure we can see that in Pb the mixing amplitudes are about 
10 times as large as those in matter with fp = 0.49 (see Figs. ^ and ^). Furthermore, 



fp 

5 ' 



in neutron stars the p-u mixing amplitude in the L mode and the a-p mixing amplitude 
are about 100 times larger than those for fp = 0.49. This fact is consistent with the result 
shown in Fig. ^. 

In summary, using QHD we have studied the p-u, a-6, a-p and 6-uj mixing amplitudes 
in isospin asymmetric nuclear matter. We have shown that the mixing amplitude is quite 
sensitive to the proton fraction and nuclear density. In particular, in neutron stars and/or 
neutron rich nuclei (like Pb) the CSB meson mixing amplitudes become considerably large 
compared with those in typical mirror nuclei. Even in light nuclei, if /„ is far from 0.5 
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(like halo nuclei pO[), the CSB meson mixing would play an important role in the nuclear 
structure. In the present calculation, we have not calculated the CSB potential generated by 
the meson mixing because, in addition to the mixing amplitude, it is necessary to know the 
CSB meson-N vertex to draw any definite conclusions on the potential [jlOl- Such a complete 



study requires more elaborate calculations and it is beyond the scope of this Letter. It would 
clearly be very interesting and important to calculate it in the future. 

We would like to thank A. Suzuki, A.W. Thomas and A.G. Williams for valuable dis- 
cussions. 
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FIG. 1. Amplitude of the p-uj mixing {pB = Po) va. the spacehke region. We take go = and 
vary \q\{= qg)- The top panel is for the L mode while the bottom one is for the T mode. The 
solid (dashed) [dot-dashed] curve represents the result with fp = 0.5(0.49) [0.51] while the result in 
vacuum is shown by the dotted curve in the top panel. 
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FIG. 2. Peak structure in the transverse mixing amplitude {pB = Po and fp = 0.51). The 
soUd (dashed) curve represents the mixing amplitude due to the proton (neutron) loop while the 
difference between proton and neutron contributions are shown by the short-dashed curve. 
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FIG. 3. Ratio of the p-ui mixing amplitude at the uj pole in asymmetric nuclear matter to that 
in vacuum. The solid (dashed) [dot-dashed] curve is for pQ (2po) [3po]- 
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FIG. 4. Same as Fig. |l]but for the a-5 (top panel), a-p (middle panel) and 6-uj (bottom panel) 
mixing amplitudes {pB = Po)- 
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FIG. 5. Meson mixing amplitudes in Pb (top panel) and neutron stars (bottom panel). We 
take Pb = 0.15 fm~^ and fp = 0.4 for Pb while pB = 3po and fp = 0.15 for neutron stars. The L 
(T) mode of the p-u mixing is shown by the solid (long-dashed) curve. The short-dashed (dotted) 
[dot-dashed] curve is for the a-S {S-lo) [a-p] mixing. 
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